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ABSTRACT 

The large-scale spatial distribution and temporal variability of cloud liquid water path (LWP) over the 
world’s oceans and the relationship of cloud LWP to temperature and the radiation budget are investigated 
using recent satellite measurements from the Special Sensor Microwave /Imager (SSM/I), the Earth Radiation 
Budget Experiment (ERBE), and the International Satellite Cloud Climatology Project (ISCCP). Observa- 
tions of cloud liquid water on a 2.5° X 2.5° grid are used over a 53-month period beginning July 1987 and 
ending in December 1991. 

The highest values of cloud liquid water (greater than 0.13 kg m -2 ) occur largely along principal routes of 
northern midlatitude storms and in areas dominated by tropical convection. The zonally averaged structure is 
distinctly trimodal, where maxima appear in the midlatitudes and near the equator. The average marine cloud 
LWP over the globe is estimated to be about 0J 13 kg nT 2 . Its highest seasonal variability is typically between 
15% and 25% of the annual mean but in certain locations can exceed 30%. Comparisons of cloud LWP to 
temperature for low clouds during JJA and DIF of 1990 show significant positive correlations at colder tem- 
peratures and negative correlations at warmer temperatures. The correlations also exhibit strong seasonal and 
regional variation. 

Coincident and collocated observations of cloud LWP from the SSM/I and albedo measurements from 
the Earth Radiation Budget Satellite (ERBS) and the NOAA-IO satellite are compared for low clouds in 
the North Pacific and North Atlantic. The observed albedo-LWP relationships correspond reasonably well 
with theory, where the average cloud effective radius (r,) is 11.1 and the standard deviation is 5.2 pm. 
The large variability in the inferred values of r e suggests that other factors may be important in the albedo- 
LWP relationships. In terms of the effect of the LWP on the net cloud forcing, the authors find that a 0.05 
kg m -2 increase in LWP (for LWP < 0.2 kg m -2 ) results in a -25 W nT 2 change in the net cloud forcing 
at a solar zenith angle of 75°. 


1. Introduction 

The importance of cloud liquid water in global cli- 
mate has been a topic of recent debate. In particular, 
the cloud-climate feedback mechanism involving 
cloud liquid water is thought to be potentially signifi- 
cant to the problem of global climate change (Paltridge 
1980). Several modeling studies have emphasized the 
temperature dependence of cloud liquid water as being 
the driving force behind this feedback (e.g., Charlock 
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1982; Somerville and Remer 1984; Betts and Harsh- 
vardhan 1987 ) . Some of these studies used in situ mea- 
surements of cloud liquid water over limited continen- 
tal regions to conclude that cloud liquid water might 
provide a relatively large negative feedback on climate. 
However, the latest modeling efforts using general cir- 
culation models (GCMs) (e.g., Taylor and Ghan 1992) 
suggest that the cloud liquid water feedback may be 
more neutral since cloud longwave radiative effects 
tend to compensate for the increased albedos. Also, 
Twomey (1991) has raised doubts about whether a 
strong dependence of liquid water content on temper- 
ature actually exists, following a reinterpretation of the 
in situ measurements. Furthermore, even if the adi- 
abatic liquid water content were to increase with tem- 
perature, there is no assurance that the clouds would 
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possess greater liquid water in a warmer atmosphere 
since other, nonadiabatic, microphysical processes 
might offset the adiabatic increase in liquid water, such 
as increased precipitation efficiency (Twomey 1991). 
Other evidence contrary to the idea that cloud liquid 
water generally increases with temperature comes from 
recent satellite measurements of cloud optical depth 
from Tselioudis et al. (1992), which demonstrate that 
optical depth and temperature are anticorrelated over 
most areas of the globe. 

Cloud feedback processes in general have been 
identified as one of several components of the climate 
system that must be understood before we can reli- 
ably predict climate change. Unfortunately, there are 
presently large gaps in our knowledge about such 
fundamental questions as, What is the global distri- 
bution of cloud liquid water and ice water? What are 
the main physical processes that determine their dis- 
tribution? and even more importantly, What are their 
respective radiative effects? Of course one of the ma- 
jor reasons for these deficiencies is that global ob- 
servations of cloud properties have only recently be- 
come available. In an effort to address these and 
other related issues, the Global Energy and Water 
cycle Experiment (GEWEX) has proposed strategies 
to monitor cloud water and other cloud properties on 
a global scale. 

A full understanding of the cloud liquid water feed- 
back and to what extent it will ultimately affect climate 
and climate change will perhaps be possible only 
through the use of GCMs. Nevertheless, observations 
can play a key role in improving our understanding of 
the factors that govern cloud liquid water and also pro- 
vide descriptions of the interactions of cloud liquid wa- 
ter with other variables of the climate system. These 
descriptions can also serve as sources of parameteriza- 
tion for simple climate models or as tests for the new 
generation of GCMs that explicitly include the predic- 
tion of cloud liquid water. 

Observations of cloud liquid water from satellite- 
based systems have been restricted almost exclu- 
sively to total column estimates using optical (e.g., 
Lin and Rossow 1994) and passive microwave meth- 
ods (e.g.. Petty 1990; Greenwald et al. 1993). De- 
spite their limitations these types of measurements 
can be extremely useful in understanding the hori- 
zontal distribution of cloud liquid water [ information 
about the detailed vertical structure of cloud liquid 
water from satellite measurements may be possible 
in the future using cloud profilers, such as the 94- 
GHz radar proposed for TRMM-2 (Tropical Rainfall 
Measuring Mission)]. The chief weaknesses of the 
optical methods include the effects of ice clouds and 
multilevel clouds and the nonlinearities arising from 
cloud brokenness /inhomogeneities within the field 
of view (FOV) of the instrument. The microwave 
techniques by comparison represent a more direct ap- 
proach but are less accurate in areas of precipitating 


cloud systems, involve relatively coarse resolution 
measurements, and are usually restricted to water 
surfaces, though retrievals may be possible over land 
when combined with IR measurements (Jones and 
Vonder Haar 1990). With the emergence of new 
spacebome microwave instruments of improved cal- 
ibration and high stability, such as the Special Sensor 
Microwave/Imager (SSM/I), it has become possible 
to monitor successfully cloud liquid water on a 
global scale since July of 1987. 

This study presents a multiyear survey of cloud liq- 
uid water over the oceans and examines the relation- 
ships of this bulk cloud property to atmospheric tem- 
perature and the Earth’s radiation budget using obser- 
vations from the SSM/I, the International Satellite 
Cloud Climatology Project (ISCCP; e.g., Rossow and 
Schiffer 1991), and the Earth Radiation Budget Ex- 
periment (ERBE; e.g., Barkstrom 1984). Our primary 
aim is to address some of the scientific concerns out- 
lined by GEWEX. These include describing the spatial 
distribution and temporal statistics of cloud liquid wa- 
ter and also providing insight into the factors that shape 
its distribution. Moreover, we wish to determine the 
coupling of cloud liquid water to the temperature field 
and the radiation budget in order to provide important 
clues as to the role of cloud liquid water in global cli- 
mate. 

A similar temperature sensitivity analysis was con- 
ducted by Tselioudis et al. (1992) using observations 
of cloud optical depth. Assuming that variations of op- 
tical depth are primarily due to changes in liquid water 
path, the results of Tselioudis et al. can be used to infer 
the relationship of cloud liquid water to temperature. 
The advantages of this study in determining this rela- 
tionship are the use of measurements of cloud liquid 
water and the elimination of possible effects from ice 
clouds at middle and high latitudes during the winter. 
In addition, this study differs from earlier comparisons 
of liquid water and the radiation budget, such as by 
Stephens and Greenwald (1991), in that we collocate 
instantaneous measurements of cloud liquid water and 
the radiation budget and focus only on low-level 
clouds. 

The first part of this paper gives a description of 
the cloud liquid water datasets and includes an esti- 
mate of their expected errors. We also explain some 
of the shortcomings of these types of datasets. In sec- 
tion 3 we describe the spatial and temporal charac- 
teristics of cloud liquid water across the world’s 
oceans over a 53-month period. Section 4 deals with 
the global, regional, and seasonal aspects of the re- 
lationship of cloud liquid water to temperature. In 
section 5 we explore the relationships of cloud liquid 
water to albedo and compare these observed rela- 
tionships to theory in selected regions. In addition, 
we investigate the effect of cloud liquid water on the 
net cloud forcing. Finally, in section 6 we present the 
important conclusions of this work. 
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2. Cloud liquid water data 

The cloud liquid water datasets are derived from 
SSM/I measurements on the Defense Meteorological 
Satellite Program (DMSP) F-8 satellite over a nearly 
continuous period from 9 July 1987 to 31 December 
1991 (the entire month of December 1987 is missing 
due to a shutdown of the instrument). The SSM/I is a 
seven-channel polarimetric radiometer that detects ra- 
diation at 19.35, 22.235, 37, and 85.5 GHz (see Hol- 
linger et al. 1990 for a detailed description of the SSM/ 
I). The datasets used in section 3 of this study consist 
of monthly mean quantities on a 2.5° X 2.5° grid. 

The physical model used to retrieve vertically inte- 
grated cloud liquid water simultaneously with water 
vapor over water surfaces is a modified form of the 
method outlined by Greenwaldet al. ( 1993). This tech- 
nique is analogous to the physical methods used in the 
retrieval of water vapor and liquid water from ground- 
based microwave radiometers (see, e.g., Westwater 
1978). The present model is also conceptually similar 
to the physical models used, for example, by Wentz 
(1983) and Petty and Katsaros ( 1992). 

An improved analytic form for the upwelling vertical 
(v) or horizontal (h) polarization brightness tempera- 
tures at the top of the atmosphere, based on a reworking 
of the results of Tjemkes et al. (1991), can be ex- 
pressed at frequency u as 

W = T s [l - ( 1 - e„ /A .„)? 2 ] - T[Z t 

“b ( i — €vih,v)Z. u tjj\ + r cos ( i — 7 2 , (i) 

where 

^ v ~ {3 oxj/ H w t w M + Hi t (2) 

*3 

= -^ [ 3 0XU H W ( t + crlJn) + H,T t JJ lM ] (3) 

are the atmospheric terms, T s is the surface skin tem- 
perature, fi is the cosine of the viewing angle, e v/lt „ is 
the surface emissivity at v- or h-polarization, % is the 
total transmission of the atmosphere CU W , V %,JJ 0XM ), r Wi „ 
is the water vapor optical depth (defined as the water 
vapor mass absorption coefficient, times the in- 
tegrated water vapor), T is the temperature lapse rate, 
H w is the water vapor scale height, and T cos is the cosmic 
background temperature. The transmission terms for 
water vapor, oxygen, and liquid water are defined, re- 
spectively, as 3 W ,„ = exp(-r^//i), V OI , v = exp (-r M ,„/ 
M) and &ij, = exp(-r^//i), where r tj/ is the liquid 
water optical depth (defined as the liquid water mass 
absorption coefficient, K tM , times the integrated liquid 
water) and t oxu is the oxygen optical depth. The con- 
tribution from the effective height of the liquid water 
is included in (2) and (3) by assuming that the cloud 
is a delta function, whereby die liquid water is concen- 
trated at a height H t . The contribution from the cosmic 


background, although very small, is included as the last 
term in ( 1 ) . Finally, the quantity cri j fi is a second- 
order term in the expansion of Z 7 that was required as 
a result of discrepancies of the analytic expression with 
more detailed calculations of T VM . The expansion yields 
a value for c of 0.25; however, a much larger value of 
3.2 was found to provide the best agreement. This value 
was used for the three lowest frequencies, but is only 
important at 22.235 GHz. 

One simplification in ( 1 ) is that the downwelling 
radiation is always specularly reflected by the sea sur- 
face, even when the surface becomes rough. This as- 
sumption is likely to be violated under windy condi- 
tions and particularly for the h-polarization component 
of the reflected sky radiation. Several investigators 
have attempted to estimate and account for the non- 
specular behavior of the reflected sky radiation (e.g., 
Wentz 1983; Petty and Katsaros 1994). This effect, 
however, is less important for the v-polarization com- 
ponent, but could be incorporated into the present 
model if needed. Another simplification is that the 
surface air temperature ( T a ) is set equal to T s . The 
expression in ( 1 ) could be easily generalized to in- 
clude T a . 

Using measurements at two different frequencies, 
the integrated water vapor (W) and cloud liquid water 
(L) can be derived iteratively from Eqs. d)-(3). 
Specification of the unknown variables and a detailed 
account of the retrieval procedure are outlined in the 
appendix. The v-polarization channels are used as op- 
posed to the h-polarization channels since they are 
much less sensitive to sea surface roughness effects for 
viewing angles typical of the SSM/I (i.e., «53°). A 
combination of measurements at 19.35v and 37v GHz 
is used for the liquid water retrievals since it results in 
lower variability in clear regions and compares better 
against surface observations than for the 22.235v and 
37v GHz pair. Measurements at 85.5 GHz would pro- 
vide greater sensitivity to changes in cloud liquid water 
for lower liquid water and water vapor amounts. We 
did not take advantage of these measurements since, 
unfortunately, both channels failed during the 53- 
month period of interest. In future work it will be 
straightforward to include measurements at 85.5 GHz 
as well in the retrieval model. 

Estimates of the retrieval errors are derived from two 
different sources. First, using a limited set of ground- 
based microwave measurements of liquid water from 
the study of Greenwald et al. ( 1993 ) , the retrievals are 
found to have an rms difference of 0.040 kg m -2 over 
San Nicolas Island and the North Sea. Second, an error 
analysis for nonprecipitating clouds, similar to the one 
presented by Greenwald et al. (1993), shows the re- 
trievals have errors of roughly 30% on average but 
which vary depending primarily on the atmospheric 
conditions and the amount of liquid water. The water 
vapor retrievals, on the other hand, have estimated er- 
rors of about 5 kg m -2 for both the 19v/37v and com- 
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bined 19v/22v/37v methods based on an extensive 
comparison with radiosonde observations by Jackson 
and Stephens ( 1 995 ) . 

The major limitation of the cloud Liquid water path 
(LWP) datasets presented in this paper, as with all mi- 
crowave-derived cloud LWP datasets, is the large un- 
certainties of the estimates in regions of rainfall. When 
the present method (based on 37 GHz) is applied to 
measurements that have significant rainfall within the 
sensor’s FOV, the retrievals will in most cases over- 
estimate the actual cloud LWP, mainly as a result of 
increased absorption/emission by raindrops. 1 Because 
of this limitation, a simple liquid water threshold is 
applied to discriminate between “precipitating” and 
“nonprecipitating” FOVs, such as done by Greenwald 
et al. ( 1993). If the retrieved liquid water exceeds 0.4 
kg m 2 , then it is flagged as possible precipitation and 
excluded from the analysis. Although we are confident 
that most retrievals suspected of being greatly contam- 
inated by precipitation have been eliminated, an inde- 
pendent method of detecting precipitation (which pos- 
sibly takes advantage of multiple frequencies ) may be 
more desirable in future studies. One must bear in mind 
that these datasets will inevitably contain low biases in 
areas of precipitation because they exclude cloud liquid 
water in precipitating cloud systems. 

When producing cloud LWP datasets for climate 
studies, it is a common and simple practice to include 
all satellite retrievals in the averaging process, whether 
or not the FOVs are clear or cloudy. As pointed out by 
Lin and Rossow (1994), this type of averaging does 
not, strictly speaking, yield the mean cloud liquid water 
since FOVs that are free of water clouds are also in- 
cluded. In their analysis, Lin and Rossow used non- 
coincident (within 1.5 h) IR data to aid in screening 
out clear FOVs. 

How one averages the LWP data depends on the ap- 
plication that is of interest to the scientist. For example, 
when comparing observations to cloud liquid water 
fields generated by GCMs, one might want to include 
clear FOVs in the average. However, if one is specifi- 
cally interested in cloud properties, then it might be 
better to exclude clear FOVs. Thus, it is highly rec- 
ommended that both types of datasets be produced. 

Unless specifically mentioned, the cloud LWP re- 
sults presented in this study include only cloudy FOVs. 
Since for our purposes it would be clearly impractical, 
given the volume of data, to use IR data to aid in iden- 
tifying clear FOVs, we use instead a liquid water 
threshold. This is similar to the scene identification ap- 
proach of Bauer and Schluessel ( 1993), who used liq- 
uid water retrievals below 0.05 kg m~ 2 as an indicator 
of clear FOVs. Here, we use a threshold of 0.048 


1 This is ignoring, of course, the separate issue of beam filling, 

which leads to underestimates of the cloud LWP. 


kg m -2 , which is three times the standard deviation of 
the retrievals from the clear sky dataset described by 
Greenwald et al. ( 1993). 

This threshold was selected for two reasons. First, it 
covers nearly the entire spread of the clear sky retriev- 
als. Second, monthly mean values of cloud amount de- 
rived from the liquid water data using this threshold are 
reasonably consistent with ISCCP cloudiness data. The 
threshold works well under most conditions, except in 
regions of rough seas and low water vapor as deter- 
mined from comparisons of the liquid water variability 
within 1° x 1° grid boxes. Thus, for W < 12 kg m~ 2 
and u > 8 m s _l , where u is the retrieved near-surface 
wind speed, the threshold is set to 0.024 kg m -2 . 

The major disadvantage in using a liquid water 
threshold is that the results will, to a certain extent, 
depend on the chosen threshold (Lin and Rossow 
1994). To test the sensitivity of the results to different 
thresholds, we compared monthly mean quantities of 
cloud LWP computed at the threshold pairs of 0.048/ 
0.024 kgm~ 2 and 0.032/0.016 kg m -2 . The latter 
threshold pair was chosen as an approximate lower 
limit since smaller thresholds yielded unrealistically 
high monthly mean cloud amounts. We find that the 
differences in the results between these sets of thresh- 
olds varied greatly depending on the cloudiness con- 
ditions. As expected, the smallest differences occurred 
in highly cloudy areas (>80% cloudiness), being less 
than 10%. In areas of 50% cloudiness the differences 
ranged from about 10% to 20%, whereas for the lowest 
cloudiness (<30%) the differences were the greatest, 
sometimes as large as 30%. 

Despite the difficulties in identifying cloudy SSM/I 
FOVs using only a liquid water threshold, the results 
presented in the following section can still provide use- 
ful information as long as it is recognized that the larg- 
est uncertainties («30%) occur generally in those 
regions with the lowest cloudiness. Obviously, there is 
a need for an independent method of detecting cloudy 
SSM/I FOVs and, even more importantly, determining 
the degree of cloudiness within the sensor’s FOV. One 
possible candidate is visible and IR measurements from 
the Operational Line Scanner (OLS) on the DMSP sat- 
ellites, which might be useful in the future processing 
of cloud liquid water datasets. 

3. Spatial distribution and temporal variation of 
cloud liquid water 

Figure 1 shows the mean SSM/I cloud LWP for the 
53 monthly averages. This figure suggests that the hor- 
izontal distribution of cloud LWP is in many ways tied 
to large-scale dynamical patterns. For example, an in- 
dication of the path of northern midlatitude storms can 
be observed, where the liquid water is often greater 
than 0.125 kg m -2 . Also, the signatures of the inter- 
tropical convergence zone (ITCZ) are evident in the 
Pacific and Atlantic, and high LWPs of over 0.125 
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CLOUD LIQUID WATER (kgm~ 2 ) 

Fig. 1. Global distribution of the mean cloud liquid water path (kg m~ 2 ) over the oceans derived from the SSM/I and 
averaged from July 1987 to December 1991 on a 2.5° x 2.5° gnd. Gray areas denote regions of land or sea ice. 


kg m 2 are common in the Indonesian area and in the 
eastern tropical Indian Ocean. 

In other regions of interest, the South Pacific con- 
vergence zone (SPCZ) is well represented by a relative 
maximum, and there is a hint of the Southern Hemi- 
sphere storm track in the southern Indian Ocean. In 



LATITUDE (deg.) 

Fig. 2. Zonal average liquid water path (kg m -2 ) averaged from 
July 1987 to December 1991. Solid curve includes only cloudy fields 
of views (FOVs) in the average, while dotted curve includes cloudy 
and clear FOVs in the average. 


areas typically dominated by subtropical highs and sub- 
sidence, the cloud LWP is generally lower except in 
those regions off the west coasts of the major conti- 
nents, particularly near North and South America, 
where extensive and persistent stratocumulus reside. 
The lowest liquid water (less than 0.09 kg m" 2 ) is as- 
sociated with areas of subsidence and generally con- 
centrated near land masses except perhaps in the central 
Pacific. 

The zonal average characteristics of the mean cloud 
LWP are shown in Fig. 2. These calculations do not 
reflect the fact that there is more ocean relative to land 
in the southern oceans than for the northern oceans. As 
a result, features in the Northern Hemisphere will be 
slightly emphasized. The zonal profile, as seen by the 
solid curve, is represented as three distinct peaks. The 
maxima in the midlatitudes are likely a result of the 
cloudiness associated with baroclinic storm systems. 
The maximum at about 5°N is caused mainly by the 
ITCZ in the eastern Pacific. Regions of minima are lo- 
cated in the subtropics. There is also an overall decline 
in the cloud LWP as one approaches the highest lati- 
tudes. The global average (weighted by area) is 0.113 
kg m' 2 . 

Also depicted in Fig. 2 for comparison is the 53- 
month mean LWP that also includes clear satellite 
FOVs. If a given FOV was determined to be clear, then 
the retrieval was set to zero to reflect the fact that no 
liquid water should be present. The differences be- 
tween these two methods of averaging is dramatic in 
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both magnitude and latitudinal shape. The “cloudy- 
clear” dataset is often nearly a factor of 2 smaller than 
the “cloudy” dataset and has four, as opposed to three, 
distinct maxima. The global average for these data is 
0.0608 kg m’ 2 . 

The seasonal character of the spatial distribution of 
cloud LWP is shown in Figs. 3 and 4. A distinct max- 
imum, which roughly coincides with the principal track 
of cyclones (see, e.g., Wallace et al. 1988), is clearly 
visible over the northern oceans during the DJF, MAM, 
and SON seasons. The large production of cloud liquid 
water mass off the east coasts of China and the United 
States, especially during DJF and MAM, is due to the 
significant fluxes of latent and sensible heat and sub- 
sequent convective instabilities that occur as a result of 
the advection of cold, dry continental air over the 
warmer ocean. 

One of the most striking features in the seasonal 
maps is the broad region of high cloud LWP ( from 0.13 
kg m" 2 to more than 0.15 kg m -2 ) during JJA in the 
North Pacific and Atlantic in Fig. 3. The conditions 
responsible for the formation of cloud liquid water dur- 
ing the summer in this region are vastly different than 
in other seasons. At this time of the year there is di- 
minished storm activity and a far greater occurrence of 
stratus clouds. The specific cause of the increased stra- 
tus cloudiness is not well understood; however, Klein 
and Hartmann (1993) have shown that the cloudiness 
is well correlated with the presence of low-level term 
perature inversions and high stability in the lower tro- 
posphere. 

To investigate further the relationship between cloud 
LWP and cloudiness, we show in Fig. 5 a plot of 
monthly mean LWP (using the cloudy-clear datasets) 
versus total cloud amount derived from ISCCP for the 
summer months from July 1987 to June 1991 over the 
North Pacific (50°-35°N; 155°-125°W). Each data 
point represents a 2.5° X 2.5° grid box average. The 
relationship between these quantities is shown to be 
reasonably well defined and nonlinear. These results 
support the idea that the increase in the liquid water is 
likely attributed, in large part, to greater cloudiness. 
However, we cannot rule out that some of the observed 
increase in cloud liquid water may also be caused by 
an enhanced thickness of the clouds. 

Another distinguishing phenomenon in the seasonal 
maps is the zonal shift of the maximum cloud LWP 
along the ITCZ. While the absolute liquid water esti- 
mates along the ITCZ are likely to be dubious, the rel- 
ative changes from one season to another may be fairly 
realistic. The most salient feature in the eastern Pacific 
is the movement of the LWP maximum to the west 
during the MAM season in Fig. 4. This result is con- 
sistent with similar shifts in rainfall maxima that were 
seen from the climatic rainfall data observed by Spen- 
cer (1993) based on measurements from the Micro- 
wave Sounding Unit (MSU) on the NOAA satellites. 


Aside from the obvious maximum of the ITCZ, other 
prominent characteristics in the Tropics can be seen for 
JJA in Fig. 3. For example, evidence of the southern 
Asian monsoon is observed as liquid water maxima 
exceeding 0. 13 kg m~ 2 on the western side of India and 
along the coast of Burma. In contrast, the cloud LWP 
is a minimum (less than 0.085 kg m~ 2 ) in this same 
region during DJF, which is attributed to the reversal 
of the monsoonal flow and a dominance of the descend- 
ing branch of the Hadley circulation. There is also an 
indication of a cloud LWP maximum in the West Pa- 
cific where the convection is known to be most vigor- 
ous. During DJF when the greatest convective activity 
advances to the southern parts of Indonesia and just 
north of Australia, the LWP values are greater than 0.13 
kg m“ 2 . We should stress that the magnitude of the 
cloud LWPs in these regions should be treated with 
some suspicion since the retrievals are less reliable 
where precipitating cloud systems prevail. 

The variation of cloud LWP with season in terms of 
zonal-mean profiles are shown in Fig. 6. All seasons 
retain, for the most part, a basic trimodal structure as 
seen earlier for the 53-month mean. In both hemi- 
spheres, particularly in the north, the equatorward push 
of cloud liquid water during the winter is observed. It 
is not clear, however, whether the seasonal differences 
in the midlatitudes are strictly a consequence of dynam- 
ical influences. The conversion of liquid water to ice 
water may also be partly responsible for the lower liq- 
uid water values during the winter season. The ITCZ 
shows little seasonal variation in its position and re- 
mains confined to the Northern Hemisphere. Evidence 
of the summertime convection in the Southern Hemi- 
sphere Tropics can also be seen as a relative maximum 
during DJF near 10°S. 

Further information about the temporal variation of 
cloud LWP can be gathered by calculating the sample 
interannual and intra-annual standard deviations, which 
for this study are defined, respectively, at each gridbox 
as 


1 !2 f i < v 

° rinler ~ 12 ^ N — 1 ^ ^ ~ Lj)" 



11/2 


(L k - L) 2 


where N is 4 or 5 years depending on the month, L, is 
the average observation over N years for a month of 
the yeary, and L is the average observation over all 53 
months. 

Calculations of <7 inter and cr intra are shown as global 
distributions in Fig. 7a and 7b, respectively. Regions 
of the highest interannual variation in Fig. 7a of over 
0.02 kg m -2 occur off the western coast of Africa near 
Cape Verde, off the northwest coast of Australia, near 
Argentina, and along the equator in the mid-Pacific. 
Other relative maxima are evident, for example, in the 
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Fig. 3. Global distributions of cloud liquid water path (kg m -2 ) over the ocean derived from the SSM/I for (a) JJA and (b) SON. 


Gulf of Mexico and in the Mediterranean Sea. Minimal 
year to year variability (less than 0.01 kg m -2 ) is seen in 
areas where high pressure systems commonly occur and 
throughout the midlatitude southern oceans. In certain lo- 
cations, such as near the northwest Australian coast and 
near Cape Verde, the interannual variations can often 
reach 25% of the annual mean. Typically, however, these 
variations are about 10% of the annual mean. 


The intra-annual variability shown in Fig. 7b, in 
many ways, resembles the interannual variability, but 
there are important differences. For instance, areas of 
high variability (as large as 0.027 kg m“ 2 ) exist in the 
tropical Atlantic, which is a result of the significant 
meridional excursions of the ITCZ. There are also large 
seasonal changes in the far eastern Pacific ITCZ. The 
Asian and Australian monsoons are exhibited as 
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Fig. 4. Same as Fig. 3 except for (a) DJF and (b) MAM. 


regions of very high seasonal swings in cloud LWP, 
sometimes as great as 30% of the annual mean. 

Another region of large (greater than 0.025 kg m -2 ) 
and widespread variability extends along the eastern 
Asian coast into the far North Pacific and Bering Sea. 
The behavior is similar just south of Greenland, though 
to a slightly lesser degree. The large seasonal variabil- 
ity in these regions may be a result of continental in- 
fluences. This is in contrast to the middle latitudes of 


the Southern Hemisphere where the impact of land 
masses is far less and the seasonal variability is smaller 
(less than 0.016 kg m~ 2 ). The lowest intra-annual vari- 
ability (less than 0.0125 kg m“ 2 ), on the other hand, 
is generally located in regions whose climate is influ- 
enced by high pressure systems. When zonally aver- 
aged, the intra-annual variability varies from 10% to 
about 20% of the annual mean. In the monsoon regions 
and parts of the midlatitudes, the variability can be 30% 
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LWP (kg m' 2 ) 


Fig. 5. Scatter diagram of monthly mean ISCCP total cloudiness 
{%) versus SSM/I liquid water path (kg m~ 2 ) for the summer months 
from July 1987 to June 1991 over the North Pacific. 


of the annual mean. The Northern Hemisphere also 
shows generally greater variability than the Southern 
Hemisphere. 

Global and hemispheric average values of marine 
cloud LWP over the 53-month period are shown in Fig. 
8 as a function of the month of the year. These quan- 
tities were derived by weighting the equal-angle grid 
boxes according to their area. The gray shading in the 
figure represents the interannual standard deviation. 
Over the globe the mean cloud LWP has a minimum 
of 0.111 kg m“ 2 in April and steadily increases to a 
maximum of 0.1 16 kg m~ 2 in December. 

When separated into hemispheric means, the 
Northern Hemisphere exhibits greater seasonal change, 
with a low value of 0.107 kg nT 2 in March and a peak 
value of about 0.1 18 kg m~ 2 in August. In the Southern 
Hemisphere, the seasonal variability is slightly less, and 
the nature of the seasonal cycle is very different than 
the Northern Hemisphere, aside from the obvious phase 
shift. For example, the lowest values of cloud LWT are 
spread out over a 4 to 5 month period and gradually 
increase to a maximum of 0.1 17 kg m~ 2 in January. 

4. Cloud liquid water temperature sensitivity 

analysis 

a. Data and analysis method 

The source of data used in the following analysis are 
the Stage Cl ISCCP cloud data products. Details con- 
cerning these datasets are described by Rossow and 
Schiffer ( 1991 ) . The comparisons are made with daily 
mean cloud-top pressure and cloud LWT data on a 2.5° 
X 2.5° grid. Low clouds are defined the same as in 


ISCCP; that is, clouds with cloud-top pressures be- 
tween 680 and 1000 mb. 

One difficulty with using daily mean gridded data is 
that although the mean cloud-top pressure may be 
greater than 680 mb, a large fraction of the grid box, 
in some instances, might be covered by higher-level 
cloudiness. As an attempt to ensure that the majority 
of clouds within the grid box are low level, we include 
only grid boxes where the low cloudiness (determined 
by ISCCP) is greater than 40%. Another potential prob- 
lem is the poor diurnal sampling of the liquid water 
observations since these measurements are derived 
from only a single polar-orbiting satellite. 

The atmospheric temperature data are also obtained 
from the ISCCP products, which are based on TIROS 
Operational Vertical Sounder (TOVS) satellite retriev- 
als. Temperatures are given at seven pressure levels: 
115, 245, 375, 500, 620, 740, and 900 mb, and at the 
surface. The daily mean liquid water data are compared 
to the layer mean atmospheric temperature averaged 
over the day, which consists of the average temperature 
from 680 mb to the surface. 

The method of analyzing the cloud LWP -tempera- 
ture relationship is very similar to that of Tselioudis et 
al. ( 1992) , so as to facilitate a comparison of the results 
of that study with the present one. The cloud LWT and 
temperature data are collected together over specified 
time periods and/or latitudinal zones and sorted into 
temperature intervals, each with a width of 15 K. The 
lower bound of the first interval occurs at -25°C, and 
each successive interval is shifted by 5 K. The temper- 
ature representing each interval is simply the average 
of the temperature data within the interval. The natural 
logarithm of the cloud LWT (L) is correlated to tem- 
perature using linear least-squares regression and the 
significance of the fit is determined from an F test. 



Fig. 6. Seasonal zonal averages of cloud liquid water path 
over the oceans for the data in Figs. 3 and 4. 
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b. Results 

Low clouds during JJA and DJF 1990 are found to 
have a globally averaged LWP of 0.095 kg m -2 . In 
terms of zonally averaged behavior (not shown), the 
largest seasonal changes in liquid water for these clouds 
occur throughout the subtropics and poleward of about 
45°N. Our analysis also suggests that wintertime low- 
level water clouds over the middle to high latitude 
oceans tend to be less optically thick than during the 
summer, while the behavior is the opposite in the sub- 
tropics. 

The relation between cloud liquid water and tem- 
perature can be characterized by the sensitivity param- 
eter defined by Somerville and Remer ( 1984) : 

1 dL = dlnL 
~ L dT~ dT ’ 

where for our application L is the vertically integrated 
liquid water content of clouds below the 680-mb level 
and T is the average temperature of the atmosphere 
below this level. 

The results of the global sensitivity analysis using 
data from JJA and DJF 1990 are given in Fig. 9. The 
quantity dlnL/dT is plotted as a function of the mean 
temperature of the 680-1000-mb layer for data ex- 
tending from 70°N to 70°S. Also shown for reference 
is the range in the sensitivity relationship one would 
expect based purely on adiabatic arguments for clouds 
of 100-mb thickness and varying height below the 700- 
mb level, where the temperature is also the layer mean. 
These calculations are based on the results of Betts and 
Harshvardhan (1987). The observations indicate that 
the sensitivities for low clouds vary greatly with tem- 
perature (indicated as circles). A quasi-adiabatic be- 
havior is observed for layer temperatures below about 
7°C, with sensitivities as high as 0.04 deg -1 . For tem- 
peratures greater than 7°C the sensitivity becomes neg- 
ative and is as low as —0.023 deg -1 for temperatures 
near 15°C. 

One source of uncertainty in the cloud LWP -tem- 
perature relationship is the impact of a different liquid 
water threshold in defining cloudy FOVs. To investi- 
gate this, the analysis was repeated using a threshold 
pair of 0.032/0.016 kg m -2 . These results are shown as 
the squares in Fig. 9. They demonstrate that the use of 
a different threshold changes the sensitivities by less 
than 0.01 deg” 1 for the global analysis. 

The relationship of d\nL/dT to temperature for low 
clouds also undergoes sizeable variations on both sea- 
sonal and regional scales. Figure 10 presents the results 
of an analysis for the same datasets but separated into 
different latitude zones of midlatitude (35° -55°), sub- 
tropics (15°-35°), and Tropics (0°-15°) for each of 
the two seasons. In the Northern Hemisphere in Fig. 
10a all regions have negative or near-zero values of 
d\nL/dT except for temperatures less than about -3°C 
(and a statistically insignificant point near 7°C) for the 


midlatitudes during DJF. The Tropics show the greatest 
sensitivity of d\nL/dT to temperature with values rang- 
ing from about -0.02 deg” 1 to -0.085 deg -1 . 

A direct comparison of the present results to those 
obtained by Tselioudis et al. (1992) is difficult since 
they used data every three hours and confined their 
analysis from 680 to 800 mb. In spite of the differences 
in these independent datasets, there exist many quan- 
titative similarities. For example, their results show that 
the sensitivities for midlatitude clouds over the ocean 
are positive (though statistically insignificant) at the 
coldest temperatures but somewhat smaller than the re- 
sults in Fig. 10a. The differences might be attributed to 
the effect of ice clouds in the observations of cloud 
optical depth. Their results also show that the sensitiv- 
ities eventually crossover into negative values at higher 
temperatures, which are as low as -0.04 kg itT 1 . These 
negative values are comparable to the values seen in 
Fig. 10a. In addition, the results of Tselioudis et al. give 
average subtropical and tropical values of about -0.05 
kg m -1 , which are generally consistent with the aver- 
age of the seasonal results in Fig. 10a, although the 
liquid water results in the subtropics tend to be of 
slightly smaller magnitude. 

The results for the Southern Hemisphere in Fig. 10b 
yield a somewhat different behavior in the temperature 
dependence and magnitudes of dlnL/dT. While the 
warmest oceans of both hemispheres are primarily as- 
sociated with negative values of dlnL/dT \ the subtrop- 
ical values are smaller in the Southern Hemisphere (on 
average about —0.01 deg -1 , but sometimes positive), 
and the tropical values are more sensitive to tempera- 
ture, ranging from 0 to almost —0.09 deg -1 . Further- 
more, the seasonal change in the southern midlatitudes 
at the coldest temperatures is far more dramatic than in 
the northern latitudes. 

Because of the coarse resolution of the measure- 
ments it is difficult to isolate the cloud LWP -temper- 
ature relationship completely from other factors that 
determine the amount of cloud liquid water over the 
grid box, such as changes in cloud thickness and cloud 
amount. The effects of variations in cloud thickness are 
reduced by using low clouds. To minimize the influ- 
ence of cloudiness on the sensitivity relationships the 
data were analyzed again but using only grid boxes 
where ISCCP total cloudiness exceeded 80%. For the 
global analysis the differences in the results obtained 
for these conditions with those presented in Fig. 9 are 
fairly small. The data become shifted systematically 
toward positive values of dlnL/dT relative to the curve 
in Fig. 9 by about 0.01 deg -1 . However, for the re- 
gional and seasonal analyses in Fig. 10, the differences 
are more significant. For example, although the basic 
temperature dependence of dlnL/dT is retained for 
most regions, the differences in magnitude of the sen- 
sitivities are sometimes as great as 0.05 deg” 1 . 

In summary, these results indicate that the LWP of 
low clouds in the Tropics and in many areas of the 
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Fig. 7. Global distributions of (a) interannual and (b) intra-annual standard deviations of SSM/I cloud liquid water path (kg m -2 ). 


subtropics and midlatitudes in the summer tends to de- 
crease with increasing temperature, while the relation- 
ship becomes generally quasi-adiabatic at colder tem- 
peratures in the midlatitudes. Tselioudis et al. ( 1992), 
in their comprehensive study, concluded that variations 
of cloud optical depth with temperature were caused 
primarily by physical processes occurring within the 
clouds. They proposed a microphysical explanation for 


the negative correlations they found between optical 
depth and temperature, ascribing them to more efficient 
precipitation processes. They also acknowledged the 
possible impact of boundary layer dynamics in strongly 
influencing these temperature sensitivities for certain 
cloud types. 

The primary objective of this section was to report 
on the observed sensitivities between cloud LWP and 
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MONTH 



MONTH 



MONTH 

Fig. 8. Areally weighted averages of cloud LWP for the globe 
(top). Northern Hemisphere (middle), and Southern Hemisphere (bot- 
tom) as a function of month. Gray shading represents interannual 
variation as one standard deviation above and below the data points. 


temperature. A thorough examination of the possible 
causes of these relationships is beyond the scope of this 
study. Clearly, the challenge in future studies will be 
to establish the relative roles of microphysical and dy- 
namical processes in determining these cloud LWP- 
temperature sensitivities. 

5. Cloud liquid water and the earth radiation 

budget 

a. Data sources and analysis method 

The top-of-the-atmosphere (TOA) broadband radi- 
ative fluxes used in the following comparisons are ob- 
tained from instantaneous measurements of the scanner 
instruments aboard the NOAA-IO satellite and the Earth 
Radiation Budget Satellite ( ERBS ). The shortwave 
(SW) channel of the scanner covers the spectral region 
from 0.2 to 5 jum and the longwave (LW) channel 


ranges from 5 to 50 //m. The ERBS scanner has foot- 
print dimensions (along-track by cross-track) at nadir 
of roughly 48 km X 32 km, while the scanner on the 
NOAA-10 is slightly larger at about 65 km x 44 km. 
The footprint size increases dramatically as the instru- 
ment looks farther away from nadir. A more complete 
description of the instrument, including its operation 
and calibration, is discussed by Kopia ( 1986). 

Collocation of the Earth Radiation Budget Experi- 
ment (ERBE) scanner and SSM/I measurements is 
achieved simply by finding those pixels whose center 
points are within 10 km of one another. Unfortunately, 
prior to 1989, the DMSP F-8 satellite experienced sig- 
nificant geolocation errors ( often more than 20-30 km) 
arising principally from ephemeris errors (Poe and 
Conway 1990). To reduce these geolocation errors we 
used the correction described by Wentz ( 1991 ), which 
is accurate to within about 10 km. 

The comparisons are confined to mostly low-level, 
summertime clouds so as to diminish the effects of ice 
particles on the albedo. The regions we have chosen to 
focus on are the North Pacific and North Atlantic. A 
number of “coincident” matches (i.e., within seconds 
to several minutes) were found in the months of June 
and July 1988 in these regions. A summary of the gen- 
eral locations, dates, and times of the matches (those 
which provided comparisons that satisfied the strict cri- 
teria described below) is shown in Table 1. 

Low-level clouds were identified by LW cloud forc- 
ing values below 45 W nrL This threshold is based on 



Fig. 9. Global relationship between d\nUdT(deg~') for low clouds 
versus the temperature (°C) of the atmospheric layer from 680 mb to 
the surface. Each symbol represents a temperature interval of width 
15 K. The observations are given for different liquid water thresholds 
used to define cloudy fields of view, i.e., the liquid water threshold 
pairs of 0.048/0.024 kg m -2 (circles) and 0.032/0.016 kg m -2 
(squares). The light shading indicates the range of the relationship 
derived from an adiabatic assumption for clouds of 100-mb thickness 
and at different heights from 700 to 1000 mb. 
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Fig. 10. Scatter diagrams of d In UdT (deg"') for low clouds as a 
function of temperature (®C) separated into midlatitude (35°-55°), 
subtropical (15°-35°), and tropical (0°- 15°) latitude regions and the 
seasons of JJA and DJF for the (a) Northern Hemisphere and (b) 
Southern Hemisphere. Filled symbols are data points with correla- 
tions below the 99% significance level. 


the results of Ackerman et al. (1992, 1994), where 
instantaneous measurements of LW cloud forcing for 
stratus clouds indicated values of up to 40-50 W m -2 . 
Clear sky LW fluxes were determined from average 
measurements, for the same month, within 1° bins as a 
function of latitude. 

To minimize the obvious effects of cloud brokenness 
on the SW measurements we chose only overcast 
FOVs, which were determined by the ERBE scene 
identification. Also, while the SSM/I’s 37-GHz chan- 
nels (alongtrack and cross-track dimensions of 37 km 
X 29 km) and the ERBE scanner have comparable res- 
olutions when the ERBE scanner looks at nadir, the 
differences in resolution can become substantial at 


large viewing angles. Thus, the comparisons are limited 
to ERBE scanner viewing angles of less than 20°. An- 
other reason for using measurements near nadir is to 
avoid the slight dependence of the albedo on the view- 
ing zenith angle, which is discussed by Suttles et al. 
(1992). Furthermore, because of the uncertainties in 
the SSM/I geolocation and in order to reduce the noise 
in the liquid water retrievals, we used 37-GHz bright- 
ness temperatures averaged from the collocated pixel 
and the four nearest pixels. 

In the comparisons, the solar zenith angles ranged 
from 70° to 86° since the DMSP F-8 satellite flies near 
the terminator. To eliminate the solar zenith angle de- 
pendence of the albedo measurements, the albedos 
were normalized to a common solar zenith angle of 75° 
using the ERBE directional model for overcast condi- 
tions. The amount of adjustment resulting from the nor- 
malization was very small for all the observations be- 
cause of the overcast conditions and since the sun an- 
gles deviated at most by 10° from the reference angle 
of 75°. The clear sky albedo, which was required for 
computing the SW cloud forcing, was determined from 
the appropriate directional model as well [see Dlho- 
polsky and Cess ( 1993) for a description of the ERBE 
directional model for clear sky ocean ] . 

The specification of the effective cloud temperature 
for the liquid water retrievals was the same as for the 
global analysis — that is, the average of the monthly 
mean surface temperature and a climatic midlevel 
cloud-top temperature. Unfortunately, coincident 11- 
fim measurements were unavailable for both the ERBS 
and NOAA-10 comparisons to aid in determining the 
cloud temperature. An attempt was also made to im- 
prove the oxygen transmission and the surface skin 
temperature in the retrieval model by using temperature 
profiles and skin temperatures derived from the TOYS 


Table I. Date, time, and approximate location of the coincident 
satellite matches in the earth radiation budget and cloud liquid water 
comparisons for 1988. 


Case 

Latitude 

ERBS and DMSP F-8 
Longitude Month 

Day 

GMT 

1 

40°-54°N 

169°E-l 66"W 

June 

1 

1719 

2 

48°-57°N 

149°- 170°W 

June 

16 

2124 

3 

54°-58°N 

139°- 154*W 

June 

18 

0540 

4 

55°-56°N 

13°-29°W 

June 

20 

2034 

5 

56°-58°N 

140°-153°W 

June 

22 

0450 

6 

55 a -58°N 

1 i°_22°W 

June 

24 

1944 

7 

40 o -45*N 

162°- 164°W 

July 

23 

1644 

8 

50°-52°N 

24°-29°W 

July 

26 

0739 


NOAA-10 and DMSP F-8 



9 

45°-58°N 

151°-157*W 

July 

7 

0512 

10 

44°-57°N 

44°-50°W 

July 

7 

2208 

11 

42 d -48°N 

29°-45°W 

July 

17 

2150 

12 

42°-58°N 

142°-149°W 

July 

18 

0434 

13 

48°-49°N 

162°-164°W 

July 

28 

0556 
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on the NOAA-9 and NOAA-IO satellites [which are part 
of the National Environmental Satellite Data and In- 
formation Service (NESDIS) TOVS sounding prod- 
ucts], Monthly mean, 2.5° x 2.5° gridded datasets of 
skin temperature and oxygen transmission were used 
in the retrievals. The oxygen transmission was com- 
puted for each sounding using the millimeter wave ab- 
sorption model of Liebe et al. ( 1993). Finally, we used 
a somewhat more conservative liquid water threshold 
of 0.35 kg m -2 to separate retrievals for nonprecipitat- 
ing and precipitating clouds. 

b. Results 

Theoretical calculations of the albedo at the TOA for 
plane-parallel clouds were done for different cloud 
LWP and for cloud droplet effective radii (r e ) of 5, 10, 
15, and 20 /xm. These values of r t are thought to rep- 
resent a typical range for most stratiform-type clouds, 
although 20 ^m is probably an extreme value. In situ 
measurements of r e in marine stratocumulus (Slingo et 
al. 1982; Stephens and Platt 1987; Albrecht et al. 1988) 
and arctic stratus clouds (Herman and Curry 1984) 
have been shown to range from about 5 to 10 /xm. 

The TOA albedo was derived from the two-stream 
model described by Stackhouse and Stephens (1991). 
The model consists of 250 SW bands from 0.28 to 3.85 
/xm and includes absorption by water vapor, C0 2 , 0 2 
and O3, and Rayleigh scattering. The 0 3 absorption is 
based on the cross-section data used in the study of 
Stamnes and Tsay ( 1990) . We implemented a modified 
discrete ordinate method to approximate the direct 
backscatter fraction as done by Stackhouse and Ste- 
phens (1991). A surface albedo of 0.21 was assumed 
for each spectral band, which was based on the ERBE 
clear sky directional model for a solar zenith angle of 
75°. The profiles of temperature and humidity were 
taken from a typical sounding over the North Pacific 
during the summer and the 0 3 density profile was ob- 
tained from McClatchey et al. ( 1972) for the midlati- 
tude summer. The cloud occupied the layer from 1 to 
3 km (889—693 mb). Lastly, the cloud optical prop- 
erties (i.e., extinction, single scatter albedo, and asym- 
metry factor) were computed for each spectral band via 
Mie theory using a representative lognormal size dis- 
tribution for marine stratocumulus. 

The results of the albedo and cloud LWP compari- 
sons for overcast, low-level clouds in the North Pacific 
and Atlantic are presented in Figs. 1 1 and 12. The com- 
parisons are shown separately for the ERBS and the 
NOAA-10 satellite because of the differences in foot- 
print dimensions of the two ERBE sensors. The obser- 
vations in Fig. 1 1 show a definite rise in the albedo 
with increasing cloud LWP for LWP less than about 
0.15-0.2 kg m~ 2 and indicate little or no sensitivity of 
the albedo to changes in LWP for larger LWPs. This 
behavior is in general agreement with theory. The re- 
lationship is similar for the NOAA-10 satellite in Fig. 

12 but is somewhat less distinct. 


Theory and observations are also shown to agree, in 
a general sense, with respect to the size of the cloud 
droplets. However, there do exist some large discrep- 
ancies in the results. A few anomalously high albedos 
(greater than 0.85) in Figs. 11 and 12 are associated 
with very large solar zenith angles (84°-86°) and in- 
dicate that the SW measurements may be suspect. 
Some of the data points corresponding to drop sizes of 
less than 5 /xm (particularly in Fig. 11) are possibly a 
result of the SW effects of ice particles, which may be 
present in some of the ERBE scanner FOVs. Also, sev- 
eral data points with drop sizes exceeding 20 /xm for 
cloud LWP greater than 0.2 kg m' 2 in Fig. 1 1 are likely 
caused by precipitation contamination of the LWP re- 
trievals. 

The mean values of r, from Figs. 11 and 12 (ex- 
cluding outlying points) are determined to be 10.6 and 
13.0 /xm, respectively, with standard deviations of 4.8 
and 6.2 /xm. Overall, the mean and standard deviation 
are 11.1 /xm and 5.2 /xm. The extent of these values are 
consistent with satellite estimates of r, for warm clouds 
in the North Pacific and North Atlantic during July 
from the study of Han et al. ( 1994), where sizes ranged 
from about 11 to 15 /xm. 

While some of the scatter in the observations is per- 
haps attributed to fluctuations in r e , the relatively large 
magnitude of these variations points to other effects 
that may be of comparable importance in influencing 
the observed albedo— LWP relationships. For example, 
uncertainties in the LWP retrievals and SW fluxes may 
cause additional variation in these relationships. Errors 
in the ERBE instantaneous fluxes (in terms of one stan- 
dard deviation) are estimated to be 12% for the SW 
and 5% for the LW (Wielicki et al. 1995) . These errors 
are dominated by uncertainties in the ERBE angular 
dependence models (ADMs), which are used to con- 
vert the measured radiances to hemispheric fluxes. The 
ADM errors may be especially important for the pres- 
ent comparison since the upwelling radiation field is 
expected to be highly anisotropic due to the very low 
sun angles. 

Macrophysical cloud effects also have the potential 
to impact the SW measurements. For instance, Cahalan 
et al. ( 1994) investigated the effects of in-cloud vari- 
ability on the albedo using fractal stratocumulus clouds 
and found plane-parallel biases of typically 10% for a 
solar zenith angle of 60°. Also, radiative transfer sim- 
ulations using a two-dimensional model suggest that 
the effects of cloud variability is most appreciable when 
the sun is low (P. Stackhouse Jr. 1995, personal com- 
munication). 

The errors associated with the two-stream radiative 
transfer model are yet another source of uncertainty in 
the albedo— LWP relationships. Two-stream methods 
are convenient because of their simplicity, but they 
have certain limitations and their performance can de- 
pend greatly on the solar geometry and the cloud op- 
tical properties. In an evaluation of several two-stream 
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Fig. il. Scatter diagram of instantaneous broadband albedo mea- 
surements from ERBS at a solar zenith angle (9 0 ) of 75° versus co- 
incident SSM/I retrievals of cloud LWP for low clouds over the North 
Pacific (filled circles) and North Atlantic (squares) during June and 
July 1988. Also shown are model simulated top-of-the-atmosphere 
albedos for different cloud drop effective radii ( r e ) based on a two- 
stream model (Stackhouse and Stephens 1991) and assuming a sur- 
face albedo of 0.21. 


methods, King and Harshvardhan ( 1986) found that for 
both conservative and nonconservative scattering the 
errors in the albedo for nearly all the methods they 
tested were within 5%-10 % for the cloud optical 
depths and solar zenith angle considered in this study. 
Finally, other less important sources of variability in 
the albedo -LWP relationships might include colloca- 
tion errors and, possibly, differences in the viewing ge- 
ometry of the SSM/I and ERBE scanner. 

Another relationship to consider, which is also of 
great significance to climate, is that between LWP and 
the net cloud forcing. The instantaneous ERBS net 
cloud forcing versus LWP is presented in Fig. 13 for 
LWPs below 0.2 kg m~ 2 . The relationship is not shown 
for larger LWPs since C nel does not vary with changing 
LWP. The net cloud forcing is defined as 

Ciel = [M)Fo( a c\T ~ Gobs) + Q.W, 

where fiQ is the cosine of the solar zenith angle, F 0 is 
the distance-corrected solar constant, a clr is the clear 
sky albedo, a obs is the observed albedo, and C LW is the 
LW cloud forcing. Also shown in this figure is the line 
determined from a least absolute deviations analysis. 
The observations reveal that the low-level clouds of the 
North Atlantic and Pacific produce a -25 W m’ 2 
change in the net cloud forcing at the TOA for every 
0.05 kg m -2 increase in their LWP. Of course, these 
results are only valid for a solar zenith angle of 75°. 
The sensitivity of C nel to changes in LWP is expected 


to be greater for the same range of LWP when the sun 
is, for example, at zenith. 

6. Conclusions 

Satellite measurements from the SSM/I and the 
ERBE scanners have been used to study the large-scale 
characteristics and radiative effects of integrated cloud 
liquid water over the oceans. The 53 months of cloud 
LWP observations exhibited many outstanding spatial 
characteristics. Some of the largest liquid water 
amounts observed (often greater than 0.14 kg m~ 2 ), 
which were possibly caused by prolific stratus cloudi- 
ness, were located in the North Pacific and Atlantic 
during the summer when these regions are under the 
controlling influence of high pressure systems and 
strong static stability. In the Northern Hemisphere dur- 
ing fall, winter, and spring the eastward progression of 
midlatitude weather systems produced a narrow belt of 
elevated cloud LWP values. Extensive and large cloud 
LWP values, often more than 0.13 kg irT 2 and some- 
times exceeding 0.15 kg m~ 2 , occurred along these 
storm pathways. In areas of the Tropics where moisture 
convergence and convection is especially strong and 
widespread the cloud LWPs were expectedly high 
(ranging from 0.13 to more than 0.15 kg m -2 ). These 
areas include the ITCZ in the Pacific and Atlantic and 
the intense convection in the western equatorial Pacific 
and Indian Ocean. The zonal distribution of cloud LWP 
consists of three primary maxima. The peaks in the 
midlatitudes are thought to be associated with storm 
activity. Another maximum in the Tropics is an indi- 
cator of tropical convection, which is caused mainly by 
the ITCZ in the Pacific. 
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Fig. 13. ERBS measurements of the net cloud forcing (W m~ l ) 
versus cloud LWP for low clouds with LWPs less than 0.2 kg nT 2 
and at a solar zenith angle (i %) of 75°. Also shown is a linear fit from 
a least absolute deviations analysis. 


The main conclusions of this study can be summa- 
rized as follows: 

(i) The spatial distribution of integrated cloud liq- 
uid water over the oceans has many characteristics- 
that are consistent with known large-scale dynamical 
features. Further research will be needed to quantify 
the effects of these dynamical processes. The average 
amount of marine cloud LWP over the globe is esti- 
mated to be 0.1 13 kg m~ 2 , while low clouds have an 
average of about 0.095 kg m~ 2 . The seasonal depar- 
tures of cloud LWP, when averaged over 2.5° X 2.5° 
regions, are typically between 15% and 25% of the 
annual mean and are largest in the Tropics and north- 
ern midlatitudes. 

(ii) The LWP of low clouds is determined to be neg- 
atively correlated with temperature throughout the 
Tropics and in parts of the subtropics and midlatitudes 
during the summer. This is in general agreement with 
previously published results based on observations of 
cloud optical depth, although there exist quantitative 
differences. At higher latitudes the correlations are re- 
versed and assume a quasi-adiabatic behavior. This re- 
sult is also consistent with the sensitivities deduced by 
Somerville and Remer (1984) based on aircraft mea- 
surements of high-latitude clouds over Russia. The re- 
lationships also exhibit significant variability on both 
seasonal and regional scales. These results cannot be 
used to imply a global cloud-climate feedback given 
that the observed sensitivities change sign depending 
on latitude. Since our results are based only on six 
months of data, many more years of data will be nec- 
essary to confirm the LWP -temperature sensitivities 
derived in this study. 


(iii) A comparison of coincident and collocated ob- 
servations of cloud LWP and albedo for summertime 
low-level clouds in the North Pacific and Atlantic re- 
vealed relationships that are, for the most part, consis- 
tent with plane-parallel theory. The mean and standard 
deviation of r e inferred from the albedo-LWP relation- 
ships are 11.1 and 5.2 [im, respectively. It is proposed 
that the large variation in r e indicates that other factors, 
including errors in the cloud LWP and SW flux retriev- 
als, may be a significant source of uncertainty in the 
albedo-LWP relationships. Further comparisons of the 
cloud LWP and the net cloud forcing for low clouds of 
LWP less than 0.2 kg m -2 show a change in the net 
cloud forcing of -25 W m" 2 for a 0.05 kg m“ 2 increase 
in LWP at a solar zenith angle of 75°. In contrast, the 
net cloud forcing for clouds with LWPs greater than 
about 0.2 kg m“ 2 exhibits no sensitivity to changes 
in LWT. 

In a follow-up study we plan to extend the albedo- 
LWP comparison to other regions and time periods. 
Future projects, such as TRMM and the Earth Observ- 
ing System (EOS), will include both radiation budget 
scanners and multichannel microwave radiometers on 
the satellite platforms. This should allow for many 
more coincident comparisons between albedo and 
cloud LWP and also provide an opportunity to deter- 
mine these relationships at smaller solar zenith angles. 
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APPENDIX 

Microwave Retrieval Method 

From Greenwald et al. ( 1993), W and L can be ex- 
pressed as 


W = 


r 1*07 ~ 19/22 


L = 


T 2 *w 19/22 T\Kwtf 


(Al) 

(A2) 


and 

<5 = **-19/22*07 ~ *w37*/ 19/22 * 

For the current study. 


r, = - Tin 


Tj — T |9v/22v + T 19/22 (W, L) 

~ AT i9v/22v 

(T s 7’cos)(l e l9v/22v) ^oxl9/22 


(A3) 
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Table A I. Model calibration factors ( K ) for the DMSP F-8, 
F-10, and F-ll satellites. Two different values are shown for A7\ 7v - 
The first values are used for the I9.35v/37v GHz retrievals; those in 
parentheses are for the 22.235v/37v GHz retrievals. 


Satellite 

AT,,, 

AT^v 

A7yrv 

F-8 

-2.2 

-2.6 

1.31 (0.95) 

F-10 

-2.2 

-2.6 

1.2 (0.84) 

F-ll 

-2.8 

-3.0 

0.93 (0.57) 



T s - Tyu - 1.5 

+ TyjjW, L) - AT 37 , 

it, - r os )(i - €, 7v )^*,„ 


(A4) 


and 


T„(W, L) = -T[Z: + (1 - OZjea (A5) 

where r, is defined for either the 19.35v or 22.235v 
channels, T i9v / 22 y and 7\ 7v are the observed vertical 
polarization brightness temperatures at 19.35 or 22.235 
and at 37 GHz respectively, l U oxl9/2 2 and & ox37 are the 
oxygen transmittances, K n9 , 22 is the liquid water mass 
absorption coefficient at 19.35 or 22.235 GHz and 
k w .i 9/22 is the water vapor mass coefficient at 19.35 or 
22.235 GHz, which represents the scaling factor be- 
tween the total water vapor optical depth and W. The 
factor of —1.5 in (A4) is the offset required for. TV,, 
computed from the simple model to agree with calcu- 
lations from a more detailed model. The quantities 
AT, 9 ,, A7V,,, and A7\ 7v account for unmodeled ef- 
fects, uncertainties in gaseous absorption, and possible 
sources of bias, including systematic biases in the 
brightness temperatures and biases in the surface emis- 
sivities. 

Values of A T 19v and A T 22 , for the SSM/I on the 
DMSP F-8 satellite are indicated in Table Al, and were 
obtained through comparisons of the water vapor re- 
trievals to radiosonde measurements. Preliminary val- 
ues for these factors for the instruments on the DMSP 
F-10 and F-ll satellites are also given in Table Al, 
which were determined by intercalibrating F-8 mea- 
surements of optical depth (r, and r 2 ) with measure- 
ments from the F-10 (coincident overpasses) and F-ll 
(overpasses within about 15-18 min). 

The physical variables in (A1)-(A5) are specified 
as follows. Here, T is set to a value of 5.8 deg km -1 
and = 2.7 K. Here H w and H, are specified as 2.0 


km and 1.5 km, respectively (note that an accurate 
value of Hi is not essential). The surface skin temper- 
ature (T,) is approximated as the monthly mean sea 
surface temperature (SST), using the datasets com- 
piled by Reynolds (1988). The water vapor coeffi- 
cients, and k w22 i are determined to be 2.3 X 10 -3 
nr kg' 1 and 6.19 x 10' 3 m 2 kg -1 . The effective cloud 
temperature, which is used to calculate x, , g , k, 22 , and 
fc, 37 (see, e.g., Greenwald et al. 1993), is estimated as 
the average of the monthly mean SST and an ISCCP- 
derived 4-year mean ( 1987-1990) midlevel cloud-top 
temperature, which was determined for each month of 
the year. This was done to improve upon the simple 
assumptions used in the study of Greenwald et al. 
(1993) and to address the concerns of Liu and Curry 
( 1993 ) and Lin and Rossow (1994). The oxygen op- 
tical depths ( t oxi/ ) depend somewhat on the tempera- 
ture profile. This dependence was parameterized 
through the surface air temperature (T a ) as T ax ,„ = c L „ 
+ Ci.JTa + c x „T\ + c 4 .„T 3 , where T a is replaced by the 
monthly mean SST in degrees Celsius. If T„ < 0°C, 
then set t ox , „ = c Uu . Table A2 shows the cubic fit co- 
efficients. 

Based on the calibration procedure described by 
Greenwald et al. (1993), k w37 was determined to be 
2.04 x 10 -3 m 2 kg -1 for the 19v/37v GHz water vapor 
retrievals and 1.93 x 10 _3 m 2 kg _l for the 22v/37v 
GHz water vapor retrievals. Values for A T 57v are 
shown in Table Al. Slightly different values of K wyi 
and AT 37v were obtained for the different methods 
since the calibration also accounts for other effects 
(i.e., unmodeled effects, surface emissivity biases, 
etc.), which may differ at 19.35 and 22.235 GHz. 

The emissivity of a smooth ocean surface was com- 
puted using the Fresnel formulas and the dielectric con- 
stants described by Klein and Swift (1977). For a 
rough ocean surface the empirical model of Petty 
( 1990) and Petty and Katsaros ( 1994) was used, where 
we excluded the emissivity bias corrections derived by 
Petty (1990). 

Retrievals of IV and L are obtained by an iterative 
process. First, initial guesses are found by setting 
T„( W, L) = 0 in (A5) and then directly computing W 
and L from ( Al ) and ( A2) . The initial guesses require 
an estimate of the surface wind speed (u), which is 
used as input into the roughness emissivity model. For 
this step, we use the wind speed retrieval algorithm 
described by Goodberlet et al. (1989). To obtain so- 
lutions for W and L, one simply uses the initial guesses 


Table A2. Cubic fit coefficients for oxygen optical depths (t,„) at 19.35, 22.235, and 37 GHz. 


Frequency 


3.86 X 10 - ’ -4.74 x 1CT 4 

4.46 X 10"' -5.57 x 10' 6 

1.33 x 1(T J -1.67 x lCT 5 


19.35 

22.235 

37 


0.0134 

0.0155 

0.0453 


6.12 x 10~‘ 
7.33 X 10~ 8 
2.26 X 10” 7 
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of W and L to compute T„( W y L) and a new wind speed 
from the following preliminary expression (G. Petty 
1993, personal communication) 

u' = u - 2.130 + 0.220IT - 4.008 x 10~*W\ 

which corrects for biases in the Goodberlet et al. al- 
gorithm that are known to exist for low and, especially, 
high water vapor. This process is repeated until the de- 
sired precision is achieved. 

For retrievals of W y the 22v/37v method gives su- 
perior results for low to moderately high water vapor. 
At the highest values of W (greater than about 55-60 
kg nT 2 ), the retrievals begin to underestimate the ac- 
tual water vapor as a result of the reduced sensitivity 
of the 22.235 GHz channel to changes in water vapor. 
Consequently, a weighted combination of the 19v/37v 
(W ]9 ) and 22v/37v GHz (W 22 ) retrievals is used. If W 22 
is greater than 53.4 kg rrT 2 , the following weighting is 
applied: 

W = W 22 + 1,75 - 0.263W 19 + 4.31 X 10~ 3 W 2 ]9 . 

This equation was determined from a least-squares fit 
of the retrievals of W Z2 and W\ 9 . When compared to 
radiosonde measurements from the study of Greenwald 
et al. (1993), the present method provides substantial 
improvements over the 19v/37v method of Greenwald 
et al. The rms differences are reduced from 1.9 to 1.4 
kg m -2 for the International Cirrus Experiment data 
and from 5.9 to 4.6 kg rrT 2 for the Tropical Cyclone 
Motion experiment data. 
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